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Abstract 

The dijet invariant mass distributions from the hadronic decay of unpolarized top quark (t — > 
bW + followed by W + — > ud) are calculated, including the next-to-leading order QCD radiative 
corrections. We treat the top decay in the complex mass scheme due to the existence of the 
intermediate state W boson. Our analytical expressions are also available in different dimensional 
regularization schemes and 75 strategies. Finally, in order to construct the jets, we use different jet 
algorithms to compare their influences on our results. The obtained dijet mass distributions from 
the top quark decay are useful to distinguish these dijets from those produced via other sources 
and to clarify the issue about the recent CDF Collaborations' Wjj anomaly. 

PACS numbers: 12.38.Bx,12.38.-t,14.65.Ha 
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I. INTRODUCTION 



Since the discovery of the top quark at the Tevatron[l[ the top quark has played a 
special role in searching for the electroweak symmetry breaking mechanism and new physics 
beyond the standard model. This can be attributed to the large mass of the top quark (about 
173 GeV ), which is almost 40 times larger than the next heaviest quark. As the Cabibbo- 
Kabayashi-Maskawa matrix element V t b approaches to 1, the top quark decays almost to 
a bottom quark and a W boson. Its decay width j^-fl is 0(GeV) } much larger than the 
typical QCD scale Aq C d ~ 300MeV, indicating that the top quark decay takes place before 
hadronization. Therefore, nonperturbative effects are not important in the properties of the 
top quark, and one can perturbatively calculate its physical quantities precisely, such as 
top quark's spin correlation. At the Large Hadron Collider (LHC) at CERN, thousands of 
top quarks are expected to be produced per year at 14 TeV. Hence a new era in top quark 
research has arrived. 

On the other hand, very recently a dijet bump around 150 GeV in the Wjj channel 
has been observed by the Collider Detector at Fermilab (CDF) at the TevatronjTJ, and it 
has attracted a 
for this anomaly 

important role in the CDF dijet excess. Moreover, the DO Collaboration reported that their 
results were consistent with the standard model's prediction in the same channel[ll|. Hence 
a careful investigation regarding the dijet in the single top production and decay is helpful. 
Even without this CDF anomaly, it is still useful to study the dijet distribution in the top 
quark decay, as a part of investigations for the top quark properties. Inspired by this, in the 
present study we will investigate the dijet mass distribution in the top quark decay. This 
work also aims at understanding the properties of top quarks. 



ot of attention. There are some explanations within the standard model 
lo| . Some studies may indicate that single top production may play an 



There are a lot of works already about top quark decays)^, 12-28|. Generally, the QCD 
eadin g or der (NLO) radiative corrections to the top quark's width amount to about 



next-to- 
-8.54% 
loop 
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17j , while the corrections of QCD two loops 2JJ, |22| and electroweak one 



23 



19[ are about —2.05% and 1.54% respectively. The nonvanishing rrib effects 
26 ] and finite width corrections l^jj reduce the Born level width by about 0.27% and 1.55%. 
All these show that the QCD NLO corrections, among others, are important for top quarks. 
Therefore, in the present study, we also put stress on the effects of QCD NLO corrections 
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to the dijet distributions in the top quark decay. 

The rest of the paper is organized as follows. Section II demonstrates the dimensional 
regularization schemes and 75 schemes. Section III tackles our scheme-independent analyti- 
cal expressions. Jet algorithms are recalled in Sec. IV and Sec. V discusses the results. The 
final section contains the conclusion. 



II. DIMENSIONAL REGULARIZATION SCHEMES AND 75 SCHEMES 

Dimensional regularization has many advantages in dealing with ultraviolet, infrared and 
mass divergences encountered in high-order calculations in a unified manner. However, 
there are still some freedoms to handle these divergences in dimensional regularization. In 
this section we recall four modern versions of frequently used schemes and adopt the first 
three in the rest of this paper. The four schemes include the conventional dimensional 



regularization 



fCDR) 



scheme (FDH) 



the 't-Hooft-Veltman scheme (HV) 29], the four dimensional helicity 
and the dimensional reduction scheme (DR)[35]. 
In CDR, only the d = 4 — 2e dimensional metric tensor is introduced, i.e. g£ = d. 
The loop momentum and the spins of vectors, regardless of whether they are "observed" 
or "unobserved" 1 are in d dimensions, whereas the spins of the spinor are in d s dimensions 
with d s > d. In this section, the observed states refer to the external states appearing in the 
hard part of the process without any subsequent hadronization. We treat d s of fermions as 
four because it is distinct from d and always appears as a global factor in computations. 



HV and FDH have many advantages in helicit y a mp 



31 



33 



itude calculations, while FDH and 



351 ]. We describe the schemes in a 



DR are two supersymmmetric preserving schemes 
unified way as explained below: 

• To maintain the gauge invariance, all momentum integrals are integrated in d dimen- 
sions. 

• The dimensions of all observed particles (hard and noncollinear external particles) are 
left in four dimensions. 



1 We call the hard and non-collinear external particles observed states and internal, soft, or collinear external 
particles unobserved states in this context. 
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• The dimensions of all unobserved particles (internal states and soft or collinear external 
states ) are treated in d s dimensions . Any explicit factors of dimension arising from 
these state should be labeled as d s temporary; these must be kept distinct from d at 
the beginning. 

We treat internal states with d > 4 in HV and FDH, whereas d < 4 in DR, i.e., all variables in 
d dimensions can be divided into a four- dimensional part and d— 4- dimensional part in CDR, 
HV and FDH, while a four- dimensional quantity can be split into d and 4 — d-dimensional 
quantities in DR. The expressions are analytic functions of d and they are continued to any 
desired regions. Setting d s = d denoted in the above items, we obtain the HV scheme, while 
setting d s = 4 results in the FDH and DR schemes. As mentioned above, the d s arising from 
dimensions of spinor space is just a global factor. Therefore, we can set this part of d s to be 
equal to 4. All of the above is summarized in Table HI 



TABLE I: Summary of dimensions in different regularization schemes. 



Regularization schemes 


CDR 


HV 


FDH 


DR 


Dimensions of momenta of observed particles 




d 






4 


4 


4 


Dimensions of momenta of unobserved particles 




d 






d 


d 


d 


Number of polarizations of observed massless vector bosons 


d 




2 




2 


2 


2 


Number of polarizations of unobserved massless vector bosons 


d 




2 


d 


-2 


2 


2 


Number of polarizations of observed massive vector bosons 


d 




1 




3 


3 


3 


Number of polarizations of unobserved massive vector bosons 


d 




1 


d 


- 1 


3 


3 


Number of polarizations of fermions 




2 






2 


2 


2 



Dimensional regularization has algebraic consistency problems with respect to 75. 75 = 
\£puvpol i *'l v l p l a , which is well defined in four dimensions. However, there are some problems 
with this definition because antisymmetric tensor e^ vp(7 lives in four dimensions only. In the 
naive definition of 75, some obviously inconsistent equalities appear. If we keep all the 
four- dimensional rules and cyclicity of the trace, the analytic continuation is forbidden 36]. 
Therefore, one should at least change one of the properties to obtain a consistent result. 
To the best of our knowledge, there are two kinds of well-known 75 strategies that have 
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been introduced: one is proposed by 't-Hooft and Veltman and proved by Breitenlohner 



and Maison 



29 



37H39| (we call it the BMHV scheme), and the other one is introduced by 



Korner, Kreimer and Schilcher 



36 



40 



41] (we call it the KKS scheme). 



As a compromise, in the BMHV scheme the anticommutation relationship between 75 
and 7 M is violated, i.e. {75, 7 M } 7^ 0. In fact, every d-dimensional quantity can be divided 
into a four- dimensional part and a — 2e part, which implies that in this scheme d > 4. 75 
anticommutes with a four- dimensional 7-matrix, while it commutes with a —2e- dimensional 
7- matrix. This definition results in some ambiguousness of chiral vector current treatment, 
e.g. 7 M i^- 7^ l^p-7 M in tree-level Feynman rules. For the current work, we take the 
symmetric version as presented in 42j, |43J], i.e. 



r 



7s 



,1+75 



— > 



1+75 



-r 



75 



2 ' 2 
1 + 75 a -75 



-r 



+ 7^- 



(1) 



2 2 ' 2 

The violation of anticommutation is also a violation of the Ward identity in axial- vector 
currents. To prevent such a violation, additional renormalization is needed [31 (Readers who 
are interested in dimensional renormalization issues can also refer to Ref s . [44M49 1 ) . This will 
be used in the next section. Although it is the first rigorously proven consistent scheme, the 
process of isolating four- dimensional and — 2e parts in the Lorentz space often suffers from 
complex practical calculations. 

On all accounts, the strategy of covariance violation in 75 has some disadvantages in 
complicated situations. On the other hand, the KKS scheme keeps the covariant anticom- 
mutations but forbids the cyclicity in the trace. In 7-matrix algebra (Clifford algebra), there 
is a unique generator, which anticommutes with all other generators in infinite dimensions. 
This generator can be defined as the 75. To avoid the cyclicity in the trace, the "reading 
point" must be chosen first, and all 75 are moved to this point before a trace is taken. This 
compromise recovers a correct anomaly as well. 

Finally, we also introduce the renormalization constants and the splitting functions in 
the CDR, HV, and FDH dimensional regularization schemes used in this paper. In order to 
avoid calculating external self energy diagrams, we choose the on-shell scheme for external 
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are 




OL s Cf 


{ 1 


4vr 


K e uv 


a s C F 




4tt 


\cuv 



5Z? S = -^( — + — -3 7£ + 31nf^-)+4 + W 



1 



5z os = _ _ (2) 

9 4vr V e c/y e IR J 

where 5Z° S , 5Z° S are on-shell(OS) wave function renormalization constants for top quark 
and light quarks, respectively, 7# is the Euler constant, and Ifdh is only nonvanishing in 



FDH scheme. The unpolarized Altarelli-Parisi splitting functions 
and CDR schemes are all listed in the following: 

P qq {z) = C F \±^-eC F (l-z), 



50, 



5l| 2 to 0(e) in HV 



Pgq{ z ) — C F e Cp z, 

z 

P gg {z) = 2N C (-?— + — + z(l - z] 
\ 1 — z z 

P q9 (z) = +[i 2 Z) -ez(l-z), (3) 



while in FDH and DR these terms should be 

1 + y' 2 

PJz) = C 



Pgq{z) = C F — . 

z 



P gg (z) = 2N C {-^- + i— - + z(l -z)^+e 2N C z(l - z), 

z 2 + (l_ z )2 

Pq 9 (z) = +[ ' Z) -ez(l-z). (4) 



III. SCHEME INDEPENDENCE AND ANALYTICAL EXPRESSIONS 



As emphasized in Sec. II, there are some degrees of freedom to regularize possible di 



vergences. Because of unitarity in QCD cross sections [51], we should expect the scheme 
independence of the well-defined physical results. In this section, analytical results are 
provided for top quark decay and subsequent hadronic decay, thus affirming the simplicity 
of these processes. Moreover, we also demonstrate that the off-shell effect in the top quark 



2 Our equations are the same as those in ref.[5l|. The discrepancies in the 0(e) parts and refs. 
were carefully discussed in ref.[51|. 
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hadronic decay is small, and narrow-width-approximation is good enough at the decay width 
level. 



A. Corrections To t —?■ bW + 

We first reproduce the well-known QCD corrections to the top quark decayj^-6|(Feynman 
diagrams generated by FEYNARTS [54[ are shown in Fig{T]). Because of the Cabibbo- 
Kabayashi-Maskawa (CKM) matrix elements 1 ~ \Vtb\ ^ \Vts\, \ Vtd\i the branching ratio 
of t — > bW + is almost 100%. For simplification, we set the CKM matrix to be diagonal 
and the mass of the b-quark equal to zero. As presented in previous works, the effect of 
nonvanishing mass of the bottom quark is negligible. Following the notations of Ref. 5j, the 
matrix element of the tree-level process t{p t ) — > b(pf,)W + (pw) with averaging over the top 
quark's spin and color is given by 



\M \ 



1 4 

e vrq 



(l -r 2 ) (l + 2r 2 ) 



(5) 



where r = and s w is the sine of Weinberg angle. We can get the leading-order width 
easily 



To 



ami 



(l_ r 2) 2 ( 1 + 2r 2) 



16s 2 „m?„ 



(6) 



where we have used the electromagnetic fine-structure coupling constant a — j-. 

To check the regularization scheme independence of these results, we first derive the 
averaged squared matrix element in the FDH regularization scheme within the naive or 
KKS 75 scheme. The virtual terms and counter-terms for renormalization are given by 

'47T/X 2 ' 



\M v \ 2 + \M ct \ 2 ) 



KKS 
FDH 



—— 2 a s Cf 
1 |-21n(l 



ln(l 



-21n 2 fl 



+ 



r 2 (l + 2r 2 ) 
2 ) - 2 Li 2 (r 2 



H_^ + 3 ln(l 
2 6 V 

ln(l -r 2 ) 



(7) 



In order to see the scheme-dependent terms, we subtract the expressions in other schemes 
by the expressions in FDH with KKS 75 treatment and use S\Ai v / c t/ re ai\ 2 = \M-v/ct/reai\ 2 — 
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{\M 



v jet I real \ j p DH 



2 )^^. These scheme-dependent terms are 



(5\M V \ 2 + 6\M ct \ 2 ) 
{6\M v \ 2 + 6\M ct \ 2 ) cl)i; 



KKS 
HV 

KKS 



Mo 12 

2 



\M 
4r 2 



4tt ' 
a s C F 



4vrr(l 



47T/i 



mi 



2\ e 



1 



1 + 2r 2 



+ 8r 2 - 1 - 8r 2 ln(l - r 2 



(5|A1„| 2 
(5|M,| 2 



<5|AU 2 ) 
5|AT ci | 2 ) 



(5\M V \ 2 + 5\M ct \ 2 ) 



BMHV 
FDH 

BMHV 
HV 

BMHV 
CDR 



lMo1 -*r> 

\Mo\ 23 -^, 
\Mo\ 2 as ° F 



47rr(l - e) V 



/ 47T/i 



2\ e 



m; 



2r 2 



4r 2 



3 + 16r 2 -8r 2 ln(l - r 2 



(8) 



These scheme-dependent terms should be canceled exactly with real corrections originated 
from soft and collinear regions. In process t(p t ) — > b(pb)W + (pw) d(Pg) > the real correction 
expressions in different schemes after integrating over the momentum of the radiative gluon 
are given by 



(\M real \ 2 ) 



KKS 
FDH 



\M \ 2 



+ 



a s Ct, 



2vrr(l 

5 
2 



2\ € 



5 - 21n(l -r 2 ) 



5tt 2 



6 



2(7r 4 - 5r 2 - 4) 
(1 + 2r 2 ) (1 -r 2 )) 



-51n(l -r 2 ) + 21n 2 (l - r 2 



2r 2 (l + r 2 )(l -2r 2 



;i + r 2 ) 2 (l + 2r 2 ) 



ln(r 2 ) + 2 Li 2 (l -r 2 



{S\M rea i\ 2 ) 

( e i a a 1 2\ J 
(0|Al rea z| ) HV 



BMHV 
FDH 



, 2 \KKS/BMHV 
I I HV 

/c-i-i-T ,o^KKS/BMHV 

(S\M real \ 2 ) CDR 



a s C F 



2"s 



0. 

4. 

| AT 12 " sCV 



47rlYl 



(5|Al u | 2 + 5|Al c t| 2 ) 

/ 4:71(1 



2\ e 



l + 2r 2 



4r 2 



-8r 2 + l + 8r 2 ln(l - r 2 



KKS 



(9) 



= -(<5|A^| 2 + 5|AU% i)(V 

Combining all the results above, we find that the results in the three-dimensional regular- 
ization schemes in the KKS 75 strategy are the same; however these are not consistent with 
the BMHV 75 scheme at present. In the BMHV 75 scheme, the violation of anticommuta- 
tion also violates the Ward identities, which is also pointed out in Ref. jsj. Furthermore, to 
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w+ 




Real 1 



Real 2 



FIG. 1: Feynman diagrams in t bW + . 



maintain the Ward identities, finite renormalization is made for axial-vector currents, 

(-pren\ ( i _ a s C p \ bare 

\ L 1*5 I FDH ~ I 2n J m5 ' 

(-pren\ | 1 _ Qs ^ F \ ybare 

\ ^JhV/CDR ~ \ tt / M5 ' 

where T m5 represents the axial- vector current. 
Thus far, we get the unique result 3 



7T 



2vr 



2tt 2 3 4 
IT ~ 2 ~ 3(1 -r 2 ) 



+ 



3d + 2r 2 



-21n( 



1 -r 2 



2x 22 - 34r 2 , . 2N 31n(l - r 2 
+2 ln(r 2 ) ln(l - r 2 ) + — — ln(r 2 ) + 



9(1 -r 2 ) 2 



l + 2r 2 



41n(r 2 ) 



+ 4 Li 2 (r 2 



9(1 + 2r 2 ) 

If we set r w 0.46, we get the well-known K factor (1 — 0.8a s 



(10) 



(11) 



3 In general, we should include finite renormalization terms of coupling constants in FDH related to con- 
ventional MS scheme [31I ] to obtain the unique physical result. However, all of our processes under 
consideration are only 0{a s ) at the QCD one-loop level. This finite renormalization is absent in our 
calculations. 
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B. Corrections To W + — > ud 



With the same procedure described in the previous subsection, we obtain the analytical 



results for the subsequent decay of W boson 



55 



56] . We labeled the momenta of the W boson, 



up (charm) quark, and down (strange) quark as Pw,Pu,Pd respectively. The diagonalization 
of CKM matrix and vanishing mass of light quarks guarantee a factor of 2 to the W boson's 
hadronic decay channel via the process W + — > ud. The diagrams of QCD correction to this 
process are all shown in Figj2j 

The lowest-order squared matrix element and decay width are 



\M 



2 e2m w 



9 ' 



4 s 2 

The contributions of virtual terms and counter-terms are 



r« = ^T- (12) 



{5\M V \ 2 + 5\M ct \ 2 ) K Hv S = ~\ M < 



{S\M V \ 2 + 5\M ct \ 2 ) B F A ^ HV = -|M 



2vrr(l - e) 








1 o ^ ' 

Z 7T 








Trr(l-e) 








1 9 nr ' 

Z 7T 









For real corrections after the phase space integration over radiative gluon momentum, we 
arrive at 



, 2 \KKS ,-rr l2 a s C F f Airfx 2 V / 2 3 17 



7T 2 



BMHV 
' FDH 



{8\Mreal\ 2 )\ 

dm l2 \KKS/BMHV |T7 l2 « s Cf l2 \KKS 
(S\Mreal\) HV = \M \ 2 — = - (S\M V f + 6\Mct \) H y • 



= -(«W + <PT*| J )£**. (14) 

After including the renormalization of the axial-vector current in the BMHV 75 scheme, 
we obtain the scheme-independent answer for the decay width of process W + — > ud, 

r = r (i + f). (15) 
10 




Real 1 Real 2 



FIG. 2: The Feynman diagrams in W + — > ud. 



C. Corrections To t — > bud 

In this subsection, we present the analytical expressions of the top quark hadronic decay. 
The corresponding graphs are shown in Figj3j As the mass of the top quark is 30 times 
larger than that of the bottom quark, we set the masses of all final states to be zero. The 
effect of the nonzero mass of the bottom quark is negligible in our results. Because of the 
intermediate-state W boson in t(p t ) —> bW + — > b(pb)u(p u )d(pd) , we treat this process in 



the complex mass scheme 53, |58( . The Born amplitude squared with averaging over the 



initial-state spin and color is given by 

w = ^^ — » > (16) 

2|s^| 4 (1 — y — z — r l y + (r z w) z 

where we have defined (p t — Pd) 2 = rn 2 y, (p t — p u ) 2 = m 2 z, r = w = ■^ L . Here, y, z, r, w 
are all dimensionless variables. We keep the width of the W boson nonvanishing. The Born 
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level decay width of this channel is 



To 



a 2 m t 
647r|s,„| 4 



[4r 2 - 2 + 2r 2 (w 2 r 2 - 3r 2 + 3) ln(? 



+r 2 (w 2 r 2 - 3r 2 + 3) ln( 



1 + w 2 



(1 — r 2 ) 2 + (wr 2 ) 2 



) 



g w 2 r 6 _ 2r 6 _ 3w 2 r 4 + g r 4 _ 1 



tan 1 (w) + tan x ( 



. wr 



71 



(17) 



By expanding it in terms of w, to 0(w°) the above result can be expressed as 



a m t 
64|sJ 4 r 2 



[(1-r 2 ) (l + 2r 2 ) w~ 1 + 



6r 4 (l - r 2 ) ) + 6r 4 - 3r 2 - 1 



7T 



(18) 



To leading order in w the result is consistent with the narrow- width-approximation (NWA), 



w 



with the Born width formulas of the top quark and W boson 



exhibited in two previous subsections. The second term is an off-shell correction, which is 
about — 0.6w relative to the first term with r 0.46. 

In naive or KKS 75 strategy within FDH regularization scheme, at QCD one loop level 4 
the squared matrix elements after renormalization with the initial-state averaged is given by 



\M V \ 2 + \M ct \ 2 ) 



2 \KKS 
FDH 



IM,| 2 



OL, Cf 



47rr(l - e) V 



6 41n((y + z)(l -y-z)) -11 



+ 



+ 6ln{{y + z){l-y-z)) 



-21n 2 (l - y-z) + 41n(l -y-z) \n(y + z) - Ahi 2 (y + z) 



2 ln(y + z) 

i-y 



A U 2 (y + z) + 7i 2 - 25 



(19) 



With the same rules as those stated in the previous subsections, the differences between 
other 75 strategies/regularization schemes and the FDH scheme in KKS 75 scheme are given 



4 Because of color flow, the W boson propagator is not involved in loops. The scalar one loop integrals with 
real masses encountered in this process were already illustrated in Ref. 



However, some analytical 

continuations should be made in calculating scalar one-loop integrals with complex arguments contrast 
to the ones with real arguments. 
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by 



(5\M V \ 2 + 5\M ct \ 2 ) 



KKS 
HV/CDR 



/ 4:71 H 

4vrr(l - e) \m\ 
3(l-y-z)(y + z)l 



|7w 12 asCp 



2\ e 



y(i-y) 



e 22/(1 -y) 



(*W + *|K*| 2 )? M " V 



(5|A4„| 2 + 5|M*| 2 ) 



' FDH 

BMHV 
HV/CDR 



(5y 2 + Tlyz + II2 2 - 5y - lis 

+4(1 - y - z){y + z) ln[(l -y-z)(y + z)})} , 

\Mo\ 2 — , 

7T 

|7Wo| 2 

4 7T 



(20) 



To check our results, we also treat the numerators of loop amplitudes in four- dimensions 
by adding the R2 terms at last. All of the results discussed above are recovered using this 
method. Because of the right-handed currents 60] of the R2 in the BMHV 75 scheme, the 
unrenormalized virtual contributions are the same within the same 75 treatment, and only 
the renormalization constants are different. 

The remaining regularization scheme- dependent terms should be canceled by the real 
radiation part. The scheme-dependent terms in real corrections resulted from the soft and 
collinear region of phase space. The two cutoff phase space slicing method given by B.Harris 
and J.Owens is used here 6JJ. The analytical result within the FDH and KKS regularization 
scheme is given by 



tm a\ KKS 
(\^\)fdb 



IM 



4tiT(1 



47T/i 5 



mi 



y z 



Q_ _ 41n[(l-y-;Q(y + z)]-ll + ^ | 

e 2 e 1 — y — z 

-2 ln 2 (l -y)-2 ln 2 (l - z) - 2 ln 2 (?/ + z) 

+4 ln[(l -y){\- z)(y + z)\ ln(f ) - 9 ln(«J c ) 

C 

-4 ln(S s ) - 12 ln(5 s ) ln{S c ) + 8 ln[(l - y - z) (y + z)} \n(8 s ) 



+61n 2 (5 s )+4Li 2 [ 



i-y 



22 - 



7tt 2 



(21) 



y z 3 

where 5 S and 5 C are two parameters to isolate the soft and collinear regions, respectively. 
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The differences between other regularization schemes and the above scheme are 



(6\M SC \ 2 ) KKS 



HV/CDR 



— 2 a s Cf 
~ |M)| A^l^e) V m t 
3(l-y-z)(y + z) l 1 

2/(1 ~y) e 2y(l-y) 
(5y 2 + 22yz + llz 2 -by- llz 

+4(1 -y-z)(y + z) ln[(l - y - z){y + z)))} 

-{5\M v \ 2 + 5\M ct \ 2 ) KKS 



HV/CDR ' 



{6\M SC \ 2 ) 
{S\M SC \ 2 ) 



BMHV 
FDH 
BMHV 
HV/CDR 



o. 

\M \ 



2 3 Oi s C p 



4 7T 



(22) 



In the BMHV 75 scheme, we also obtain the scheme-independent results after including 
the finite renormalization to the axial-vector currents. This was done in order to maintain 
the Ward identities as already shown in the last two subsections. 

In the hard noncollinear phase space region, we treat the squared matrix element of 
t(pt) ~~ >* KPb) u (Pu)d(pd)g(pg) in four dimensions. Dimensionless variables are redefined as 
follows: 



ipt ~ Pgf = m 2 x, (p t - pu) 2 = m 2 y, (p t - p d f = m\ z, 



{Pu + Pdf = m 2 t k, (p u + p g ) 2 = m 2 l,r 



m w 



-i w 

m w 



(23) 



The averaged squared amplitude is 

„4 „2 n_ ( 1 



m ,2 _ 3 e 4 gj C F 



\s w \ 4 m 2 \ (k — r 2 ) 2 + (wr 2 ) 2 (1 — x) 2 (l — y — z — k) 
[(x - 3) (k + I) 2 k + 2 (xz - 2x - y - 4z + 4) k 2 

+ ((2xz - 5x - Ay - IO2 + 11) I - x 2 - y 2 - 7z 2 + Ay - Ayz + 15z 

+ (y 2 - 2y + z 2 - 7z + 6) x - 7) k + (1 - z) (2 - x - y - z) (1 - x - 2z) 

- (x + 2y + 2z - 3) I 2 - (x + 2z - 2) {x + 2y + 2z - 3) Z] 
1 1 
(2 — x — y — z — r 2 ) 2 + (wr 2 ) 2 l(x + y + z + k + l — 2) 
[(2 - x - y - z) (1 - 2 - k - I) 2 - (k + I) (1 - z - k - I) 

+ (l-y)(l- y (2-x-y-z))]} (24) 

There are two kinds of Breit-Wigner distributions of the W boson in Eq. (p4|) . The first 
term originated from the first two real diagrams, while the second is contributed by the last 
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Real 2 Real 3 Real 4 

FIG. 3: Feynman diagrams in t — > bud. 

two final state radiative diagrams. Because of color flow, there is no interference observed 
between the first two and the last two diagrams. 



IV. JET ALGORITHMS AND PHASE SPACE 



At high energy colliders, it was pointed out that the observed jets provided a view of 



parton (e.g. gluon and quark) interactions occurring at short distances 62]. At leading- 
order (LO) level, partons can be naively treated as jets, while at NLO level this coarse 
treatment often suffers from soft and collinear divergences. Therefore, an infrared-collinear 
safe jet definition is necessary in investigating strong interaction physics. Nowadays, these 
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jet definitions play important roles in collider physics. Following the jet definition description 



in Refs. 



63 



64j, the requirements implemented in a jet algorithms are as follows: 



• simple to use in experiments and theoretical calculations, 

• infrared and collinear safe, 

• small hadronization corrections. 

At hadron colliders, a well-defined jet algorithm must be able to factorize initial-state 
collinear singularities; they should also be isolated from the contamination of hadron rem- 
nants and underlying soft events. 

Since the advent of jet production in electron-positron and hadron colliders, it has be- 
come one main tool in QCD research. Many kinds of algorithms have been proposed 



and developed. Essentially, the two classes of jet algorit 



tering algorithms 



65 



66] and the cone-type algorithms 



62 



ims present mainly the clus- 



67H69|. In the present study, 



we focused on t 
algorithm 



63 



re three popular inclusive clustering algorithms, namely the /c^-clustering 



64 



70|,the Cambridge/Aachen clustering algorithm(CA) 



71 



72 1 and the anti- 



k± clustering algorithm 
described uniformly: 



73[ respectively. These three inclusive clustering algorithms can be 



Define a distance dij = minip^j^fj)^^- between each pair of protojets i and j, as 
well as a distance dis = Pri between each protojet i and the beam,with r = +1, 0, — 1 
corresponding to k^,CA, and anti-/cj_ respectively. 

Find the smallest of all the d^ and diB and label it as d m i n . 

If d min is a dy, then cluster protojets i and j as a new protojet with a selected com- 
bination procedure. If the distance between protojet i and the beam is the shortest, 
set the protojet i aside and leave it without any further clustering as a possible jet 
candidate. 



Repeat the items above until there is no protojet left. 



Perform some cuts (as in the experiment) to select jet(s) of interest. 



Here ARij = \J (r/i — r]j) 2 + ((pi — (pj) 2 (77 and (p are rapidity and azimuthal angle respec- 
tively). As E can be measured at e + e~ colliders rather than only at hadronic colliders, 
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one should use E instead of pr and ASij = y A6f- + ^sin e ^ 3 Afa^j instead of AR^ at 
e + e~ colliders. 



It was also emphasized in Ref. [70] that traditional cone- type jet algorithms were related 
to clustering algorithms by the approximation R c i us ter = 1.35 x R cone . 

In the present study, we only used the three clustering algorithms with the E-scheme 
recombination to reconstruct our leading two jets from top quark hadronic decay in the 
next section (one can also use other recombination procedures as suggested in Ref. 63j and 
references therein). In addition, we used hadron collider clustering algorithms and electron- 
positron collider clustering algorithms but without any cut in our calculation. 

The last topic of this section is about a phase-space integration treatment. Given that we 
should reconstruct the four momenta of all final states in order to reconstruct two leading 
jets, we built up the n-particle phase space iteratively by nested integration over invariant 



masses and solid angles of outgoing particles, similar to the strategy in Ref. 



74|. 



V. RESULTS 



The dijet invariant mass distributions with different clustering jet algorithms are pre- 
sented in this section. 

As discussed in the previous section, we used two variations of clustering jet algorithms 
in our top decay process in the cm. frame of the top quark. In these two variations, 
we chose the distances defined at hadron colliders ( i.e. use px ) and e + e~ colliders( i.e. 
use E), respectively, to reconstruct the final jets. Afterward, two leading jets in energy E 
were chosen to construct their invariant mass rrijj. Here,we call the first types KT1, CAl, 
anti-KTl, while the second types are denoted as KT2, CA2, anti-KT2. 

The following input parameters are used: 

a' 1 = 129, a s {m z ) = 0.119, 

m w = 80.399 GeV,T w = 2.085 GeV,m z = 91.1876 GeV, T z = 2.4952 GeV, 

4 = 1- m \ - % mW l W = 0.222657 - 1.11098 x 10~ 3 *, (25) 
m% - % m z T z 

with two groups of top quark mass and renormalization scale fi choices, i.e. m t = 175 GeV, 
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fi = 80.4 GeV and m t = 172.5 GeV, fi = m t = 172.5 Gel/ 5 . 

We varied the parameter R in the CAl jet algorithm and compared its influence on our 
results in FigJH Only when R > 1.0 in this algorithm, the infrared and collinear safety in 
each bin is maintained. Therefore, to ensure reliability, we choose R = 1.4 in the first type 
jet algorithms and R = 1.3 in the second ones for the rest of the paper. There are some 
interesting characters in these two figures. The variation of R slightly changed our domain 
region (110-150 GeV) both in LO and NLO level. The larger R reconstructs a smaller 
number of final jets; it makes the number of events in the last bin (170-175 GeV) larger with 
larger R. At LO, the distributions dropped sharply below 110 GeV and vanish below 100 
GeV, as shown on the upper panel of FigJH In contrast, a NLO QCD correction resulted in 
the smooth descent of the low energy tail. The peak in FigJH (lower panel) between 80 GeV 
to 85 GeV is the W boson's resonance. 

Histograms in Figs]5] and |6] establish the influences of clustering jet algorithms to the 
dijet invariant mass distribution. The LO distributions reconstructed by various algorithms 
are almost indistinguishable. In comparison,there are some differences in the substructures 
of NLO histograms. The combination sequence of protojets is responsible for these tiny 
distinctions 6 . Soft protojets may be clustered before the hard ones in k±, while the situation 
may be totally different in anti-/cj_. For comparison,we also plot the histograms with m t = 
\72.hGeV and /i = 172.5GW in Figd 



5 As shown in Sec. Ill, the only scale /i dependence in 0(a 2 a s ) is a s (fi), which is just a global factor and does 
not change our dijet invariant mass distribution significantly. However, the top quark mass dependence 
in our results is much more complicated. Therefore, we choose two top quark mass benchmark points to 
investigate its influence on our curves' shape, and do not plot the scale dependence in this paper. 

6 Statistical uncertainties are also responsible for these differences in the histograms. They change our 
results by about 4 percent. 
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FIG. 4: The LO (upper panel) and NLO (lower panel) dijet invariant mass distribution from top 
decay with different R using the CA1 clustering jet algorithm (m t = 175 GeV, /i = 80.4 GeV). 
Plotted are R=1.0 (solid line), 1.2 (short-dashed line), 1.3 (dotted line), 1.4 (long-dashed line), and 
1.5 (dot-dashed line), respectively. 



19 



+ 

— : 

T 

c 



0.025 



0.020 



0.015- 



0.010- 



0.005- 



O.OOOi 



60 



0.020 



$ 0.015 



0.010 



0.005- 



O.OOOi 



60 



LOR KT1 = 1.4 
LOR C ai = 1.4 
LOR AntiKT1 = 1.4 



80 100 120 

MjjiGeV) 



140 



160 



NLOR K ti = 1.4 
NLORcai = 1-4 
NLO R AntiKT1 = 1.4 

r - 




80 100 120 

MjjiGeV) 



140 



160 



FIG. 5: The influence on distribution with different clustering jet algorithms of the first type 
{mt = 175 GeV, [i = 80.4 GeV). LO is in the upper panel, while the lower panel is for the NLO 
results. Plots are KT1 (solid line), CAl (dashed line), and anti-KTl (dotted line). 

VI. CONCLUSIONS 



We have performed QCD radiative corrections to the dijet production in the unpolarized 
top quark hadronic decay in the complex mass scheme. We carefully checked the indepen- 
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dence of dimensional regularization schemes and 75 strategies in our analytical formalism. 
Applying different clustering jet definitions, we obtained our final dijet invariant mass dis- 
tributions. The obtained dijet mass distributions from the top quark decay are useful to 
understand the top quark properties and also to distinguish these dijets from those produced 
via other sources. Therefore, these results are useful in investigating the recent CDF Wjj 
anomaly and clarifying this interesting issue. Furthermore, a more careful investigation for 
top and W boson associated production at hadron colliders will be definitely needed. 
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FIG. 6: The influence on distribution with different clustering jet algorithms of the second type 
{mt = 175 GeV, = 80.4 GeV). LO is in the upper panel, while the lower panel is for the NLO 
results. Plots are KT2 (solid line), CA2 (dashed line), and anti-KT2 (dotted line). 
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FIG. 7: The influence on distribution with different clustering jet algorithms of the second type 
{mt = 172.5 GeV, ji = 172.5 GeV). LO is in the upper panel, while the lower panel is for the NLO 
results. Plots are similar to Figj6j 
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